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Harmonic oscillations of the walls of a turbulent plane channel flow are studied by 
direct numerical simulations to improve our understanding of the physical mechanism 

■ for skin-friction drag reduction. The simulations are carried out at constant pressure 
gradient in order to define an unambiguous inner scaling; in this case, drag reduction 
manifests itself as an increase of mass flow rate. Energy and enstrophy balances, carried 
out to emphasize the role of the oscillating spanwise shear layer, show that the viscous 
dissipations of the mean flow and of the turbulent fluctuations increase with the mass flow 
rate, and the relative importance of the latter decreases. We then focus on the turbulent 
enstrophy: through an analysis of the temporal evolution from the beginning of the wall 
motion, the dominant, oscillation-related term in the turbulent enstrophy is shown to 
cause the turbulent dissipation to be enhanced in absolute terms, before the slow drift 
towards the new quasi-equilibrium condition. This mechanism is found to be responsible 
for the increase in mass flow rate. We finally show that the time-average volume integral 
of the dominant term relates linearly to the drag reduction. 



1. Introduction 

The reduction of skin-friction drag in wall-bounded turbulent flows is an important and 
challenging area of fluid mechanics. Its difficulty lies both in the extreme complexity of the 
physics underlying turbulence and in the resistance of such flows to change favourably 
when disturbed by external agents. The interest in the subject is steadily growing as 
the viscous action exerted by turbulence causes dramatic energy losses in flow systems 
of technological relevance, such as oil and gas pipelines, high-speed aircraft wings, jet 
engines intakes, and turbine blades. Even a small reduction of turbulence activity, and 
thus of wall friction, translates into improved system efficiency and therefore into lower 
fuel consumption. Further potential advantages are the attenuation of noise, structural 
vibrations, and aerodynamic heating. 



f Present address: ONERA, Departement d'Aerodynamique Fondamentale et Experimentale, 
8, rue des Vertugadins, 92190 Meudon, France. 
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Active turbulent drag reduction techniques, for which energy is introduced into the 
system, have received widespread attention, owing to the (so far) hmited performances 
achieved by most passive techniques. Closed-loop feedback control strategies represent 
an emerging field of research whe re the activation is usually applied at the wall as wall- 
normal distributed transpiration (IKim fc Bewlevll2007: Kasagi et aLll20096 ). Open- loop 
techniques, for which the control law is predetermined, usually operate at much larger 
spatio-temporal scales and do not require distributed sensing. As for this type of forc- 
ing, near-wall flows have been excited by Lorcntz forces or Dielectric Barrier Discharge 
plasma actuators, alternating wall suction and blowing, unsteady cross- flow pressure gra- 
dients and different types of wall motion with the intent to disrupt the self-sustaining 
turbulence production mechanisms. With regard to the space-time distribution of forc- 
ing, both spanwise- and streamwise-traveling waves have been employed. A recent volume 
([Leschziner et aLll201lh contains several contributions to the subject. 

We consider here the simplest amon gst such open-loop techniques, i.e. the harmonic 
spanwise wall oscillations introduced bv lJung et all 1 19921 ). The oscillating wall has been 
chosen because it can be regarded as a paradigm for a larger class of drag-reduction 
techniques and because it offers the largest amount of available experimental and nu- 
merical data and the smallest number of forcing parameters. This flow has been studied 
mainly through turbulence statistics, flow visualizations of the near-wall modified flow, 
and simplified models which attempt to explain the physics behind drag reduction. Vari- 
ous mechanism for drag reduction have been pro posed, such as the relative displacement 
of near-wall structures (Baron & Quadrio"l996") and the creation of negative spanwise 
vorticity during the oscillation cycle (Choi et al. 1998,). In spite of such efforts, the answer 
to fundamental questions, such as of why the turbulent kinetic energy and the friction 
drag decrease and how the wall forcing can be modified most efficiently to achieve the 
largest net energy saving, still remain elusive. 

The objective of the present work is therefore to gain further insight into the physics of 
an incompressible channel flow with spanwise wall oscillations. Although our conclusions 
will be limited to the oscillating wall, it is our hope that they will bear some general- 
ity. The focus is on how the energy transfer between the mean flow and the turbulent 
fluctuations is affected by the wall motion and on the role played by the forcing on the 
modification of the turbulent enstrophy. The approach is to identify those terms in the 
equations which are directly affected by the spanwise forcing and to single out the domi- 
nant one(s). Another important point is to study the energy transfer during the temporal 
evolution from the start-up of the wall motion with the aim of explaining the decrease 
of skin- friction coefficient. As statistical and flow visualization studies on drag-reducing 
flows are affected by the reference quantities used for dealing with dimensionless quan- 
tities, an important and critical choice made in this study is to carry out the numerical 
simulations with a constant streamwise pressure gradient: this provides us with a clear 
and unequivocal inner-units scaling. 

The flow conflguration, the numerical procedures, the flow field decompositions and 
the basic flow statistics are presented in SJ^J The analysis of the energy budget is given in 
331 The turbulent enstrophy budget, chosen as a convenient substitute for the turbulent 
dissipation budget, is discussed in 21 Section Sj5] contains a summary of the results. 



2. Flow configuration and numerical procedures 

An incompressible fully-developed turbulent channel flow between two indefinite par- 
allel flat plates at a distance 2h*, driven by a constant streamwise pressure gradient 11*, 
is studied by direct numerical simulations (DNS). Dimensional quantities are indicated 
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by the symbol *. The coordinates a;*, y* , z* indicate the streamwise, wall-normal and 
spanwise direction, respectively. The two walls at y* = and y* = 2h* oscillate in phase 
along z* according to w*(x* ,z* ,t*) — A* cos (27rt*/T*) , where t* is time and T* is the 
oscillation period. Quantities are scaled by viscous units, i.e. by the kinematic viscos- 
ity of the fluid v* and the friction velocity u* — \/t*^/ p*, where is the time- and 
space-averaged wall-shear stress and p* is the density of the fluid. The friction velocity 
Reynolds number is Rtr = u%h* /v* = 200. As 11* is constant, the momentum balance 
at the walls shows that, once the oscillating-wall regime is established, (and therefore 
Rcr) retains the fixed-wall value. It follows that a unique wall-unit scaling is defined. 
Since this is the only scaling used throughout the paper, we omit the customary symbol 
-|- marking inner-scaled quantities. (We would like to point out here that the advantage 
in carrying out DNS at constant pressure gradient is not general. There is a potential 
impact on the computing costs for cases, like the present one, where an abrupt change 
of one parameter is introduced. As a general rule, a constant flow rate allows the wall 
friction to reach the new state sooner, thus yielding a shorter transient. However, this 
advantage is compensated by the integration time required to obtain a reliable value 
of the mean wall friction, which is averaged over two spatial directions only at each 
time step. This quantity thus presents larger temporal fluctuations than the flow rate, 
which is a volume-averaged quantity. At the present value of the Reynolds number, the 
computational cost of the two approache s is comparable.) 

Details on the DNS code are found in Luchini fc Quadriol (|2006t ). The computational 



domain has dimensions of — 6T:h*, L* = 2h* , = Stt/i* in the three directions. 
The wall-normal direction is discretized by 160 mesh points and 320 x 320 Fourier modes 
are used along the homogeneous x* and z* directions. The time step is At = 0.1 which 
guarantees that the CFL condition is amply verified for the chosen time integration 
scheme (a three-substep low-storage Runge-Kutta). The mean velocity profile and the 
variance of velocity fluctuatio ns for the fixed-wall case, sh own in figures [T] and El have 
been compared with those bv ldel Alamo &: Jimened ( 2003f l at a slightly lower Reynolds 



number and excellent agreement has been found. When the wall oscillates, the amplitude 
(maximum speed) of the wa U motion is A = 12. The effect of A on drag reduction has 



been previously studied by iQuadrio fc Riccd (|2004l ) and is not considered here. The 



calculations span the range < T < 175, and most of the paper discusses one case with 
T = 100, which is the oscillation period unless otherwise indicated. 

2.1. Averaging operators and flow field decomposition 

This paper employs different types of space and time averages and the relevant operators 
are presented here. A quantity f(x,y,z,t) is averaged along the homogeneous x and z 
directions as 

fiy,t)^-n- / fix,y,z,t)dzdx. 

J^xJ^z Jo Jo 

The velocity and the vorticity fields, U = U(x, y, z, t) and O, — Jl(a;, y, z, t), are decom- 
posed as follows 

U= {c/(y,i),0,W?(y,t)} + {u,«,u;}, O = {f},(y, t), 0, f2,(2/, t)} + {c^„ c^^, wj, 

where (ix = {ll2)dW /dy and (l^ — ~~{l/2)dU/dy since V — 0. At statistically steady 
state, a quantity f{y,t) is averaged over the N periods of oscillation as follows 

n=0 
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Note that these averaged quantities henceforth depend on the 'window' phase-average 
time T, < T < T. Alternatively, they can be observed as a function of the oscillation 
phase (j> = 27rT/T. A quantity f{y, t) is then averaged over T according to 



/(y,T)dr. 



A global quantity [f]^ is obtained by integrating (/) (y) along y, as follows 

rh 



[/] 



(/) (y)dy. 



A transport equation is defined as global when its terms are global. All statistical samples 
are doubled by averaging over the two channel halves, by properly accounting for the 
existing symmetries. 

2.2. Definition of turbulent drag reduction 
The skin- friction coefficient is defined as C/ = 2rj^/p*C/^^, where is the bulk velocity, 



u* 



(2.1) 



Following lKasagi et al\ (|2009Qf ). the drag reduction R is defined as the change of Cf with 
respect to the fixed- wall value C/^o, i-e. R ~ (Cf.o ~ Cf)/Cffl- When 11* is constant, R 
is due to the increase of mass flow rate. As 



C 



f 



R may be written as 



R 



r/2 



^1 



(2.2) 



(2.3) 



2.3. Basic flow statistics 

Figure[T](left) shows that the mean velocity profile (u^ increases significantly throughout 

the channel for T — 100 (i? = 0.31), while the wall-shear stress remains constant, in 
agreement with e xperimental studi e s where the drag-r educed friction velocity was used 



for inner scaling (jChoi et a/.lll998l: iRicco fc Wull2004j) . Figured] (right) shows that R 
increases sharply with T up to the optimum Topt ~ 70 and then decays at a slower rate. 
This behaviour is well documented by previous numerical studies, although quantitative 
differences can be ascribed to different scaling procedures. For example, the optimum 
period at constant A* is typically reported to be Topt ~ 100 — 125 at constant mass flow 
rate if the flxed-wall u* is used for scaling. 

The variance of the turbulent velocity fluctuations and the Reynolds stress component 
(uv) are shown in figure [2] (left). The wall motion primarily affects ^'^^^ up to y « 30; 
the peak decreases and its position shifts upward from y « 14 to y « 20. The profile 



of 



is largely unvaried, while that of 



increases up to y « 40. As discussed by 

iQuadrio fc Riccol (|201l[) in the context of streamwise-travelling waves, it appears that the 
large reductions of turbulence fluctuations for all the velocity components often reported 
in the literature are largely a byproduct of the outer scaling employed to compare flows 
that in fact have different values of Rer owing to drag reduction. The Reynolds stresses 
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Figure 2. Left: Wall-normal profiles of the variance of velocity fluctuations and of Reynolds 
stress (uv'^ for fixed-wall (thin lines) and oscillating-wall (thick lines) cases. Squares: u, dia- 
monds: V, circles: w, triangles: (uv'^. Right: Wall-normal profiles of Reynolds stresses vw at 
different values of the phase angle 0. 



(uv) are attenuated up to y « 60. This is consistent with iMarusic et al\ (|2007f )'s finding 
on the relati on between drag reduc tion and a weighted integral of (uv), an extension of 
the result bv lFukagata et al. ( 2002 ) to the E-constant case. 

The wall oscillation induces the additional phase- varying Reynolds stresses vw, shown 
in figured] (right). This term is null in the fixed-wall case. At opposite phases of the cycle, 
the vw profiles show the same behaviour with opposite sign, which leads to {vw) — 0. 



3. Energy balance 

As follows from (|2.2p . the reduction of C/ can be understood by studying how Ub 
increases. As a first step, we therefore study the transport equations for the mean kinetic 

energy (MKE),(^C/2 + W^^ /2, where Ub appears explicitly, and for the turbulent kinetic 

energy (TKE), g^/2, where — UiUi. (The Einstein summation convention of repeated 
indices is adopted henceforth and the subscripts i = 1,2,3 denote the x,y,z directions 
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4> = o 



Figure 3. Wall-normal profiles of —vwdW/dy at different phases (j> (solid lines) and their 

time-averaged value (dashed line). 



and the corresponding velocity and vorticity components.) These two equations are then 
summed to obtain the global balance for the total kinetic energy. 

3.1. Mean kinetic energy balance 
The transport equation for MKE reads 

^d(u' + wA ^ Qfj d(B,u) dii^w) 

^ ^ ' C/n =m—+€S}- V / \ / 



2 dr dy dy 






/ dW 




\ dy 


9 


10 



(3.1) 



Term 1 denotes the temporal change of MKE and term 2 is the work per unit time 
done by H, i.e. the power used to drive the flow along the x direction. Thanks to the wall 
oscillation, the system absorbs more kinetic energy than in the fixed- wall case through the 
increment of U. Term 3 is the work of deformation carried out by the Reynolds stresses 
uu, through which energy is exchanged between the mean flow and the fluctuating flow. 
Term 4 indicates the work of deformation done by the Reynolds stresses ifw; similarly 
to term 3, it transfers energy between the mean flow and the fluctuating flow. Terms 3 
and 4 appear with opposite sign in the TKE equation, as shown in i j3.2l The transport 
works performed by the Reynolds stresses uv and vw are described by terms 5 and 6, 
respectively. Terms 7 and 8 are the transport works done by the mean streamwise and 
spanwise viscous stresses, respectively. Term 9 is the viscous dissipation of MKE by the 
wall-normal^radient of U, while term 10 is the viscous dissipation by the wall-normal 
gradient of W. 

The second part of term 1 and terms 4, 6, 8, 10 are directly related to the wall oscilla- 
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tion, since W appears explicitly in their expressions. The turbulent^roduction term 4, 
—vwdW/dy, which is absent in the fixed- wall case because vw and W are null, is shown 
in figure [3] at different phases of the cycle. Although it is negative during part of the 
cycle (mainly for y < 15, when it instantaneously extracts energy from the turbulent 
fluctuations to enhance MKE), it is positive for most of the cycle, i.e. its average con- 
tribution is to transfer MKE to the turbulent fluctuations; see dashed line in figure [31 

which represents — (vwdW / dy'^ . 

As the primary interest resides in the change of Ub, the MKE equation (|3.ip is now 
time-averaged and integrated along y to make Ui, appear in the energy balance. Time 
averaging eliminates term 1 because of time periodicity. Terms 2, 3, 4, 9, 10 are retained 
as is term 8 since W is non-zero at y = 0. Terms 5 and 6 disappear because uv and 
vw are null at y — and at y = h. Term 7 becomes null because U = Q at y = Q and 
dU/dy = at y^h. 

The global transport equation for MKE is 



UbT^ + ( A 



dW 
dy 




where = (^dU/dy\y^oj. The first term on l.h.s. comes from term 2 in (|3.ip and 
represents the global energy per unit time pumped into the system through the external 
pressure gradient 11. Term is the energy input given by the wall motion, and denotes 
the energy spent to move the walls against the frictional resistance of the fluid. It stems 
from the transport term 8 in p.ip . Terms Vuv and Vvw, which originate from terms 3 
and 4 in (13. ip . are a sink for MKE and appear in the global TKE balance as production 
terms. Terms Vjj and T>w, which stem from terms 9 and 10 in (|3.ip respectively, denote 
the global viscous dissipation due to the gradients of the mean streamwise and spanwise 
velocity components. Equation (|3.2p represents the first step toward understanding drag 
reduction because Ub now appears explicitly. It states that part of the energy input, 
UbTw + Sw, is transferred to the turbulence via Vuv and Vvw^ and the remaining part is 
dissipated into heat through Vu and Vw- 

3.2. Turbulent kinetic energy balance 
The transport equation for TKE reads 




d {vp) 
dy 



dy 



uv 



dy 



■ vw- 



dW 

dy 



lav 

2 9y2 



dui dui 
dxj dxj 



(3.3) 



where p is the turbulent pressure. The temporal change of TKE is expressed by term 
1, while terms 2 represent the work of transport done by the total dynamic pressure 
of turbulence. Terms 3 and 4 denote production of TKE and also appear in the MKE 
equation (|3.ip with opposite sign. Terms 5 and 6 together represent the combined effect 
of the work done by the viscous shear stresses of the turbulent motion and of the viscous 
dissipation of TKE into heat. Term 6 is often referred to as the pseudo-dissipation (see 
Popd (l2000t ) at page 132). The turbulent production term 4, vwdW /dy, is the only one 
containing W explicitly. 
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Analogously to the analysis of the MKE equation, time averaging and integration along 
y lead to the following simplifications. Term 1 disappears because of time periodicity. 
Terms 2 become null upon y-integration because of the no-slip condition a,t y — and 
vp = ~ at y = h. Term 5 is also null because 



dy \dy 




dy ^ 



0, 



as dq^ /dy ^ at y ~ h and q = at = 0. 
The global transport equation for TKE is 









uv — 


+ 


vw—— 


dy 




dy 




9 






Pi,™ 



+ 



dui dui 
dxj dxj 



(3.4) 



where Vuv and are as in (|3.2p . The next-to-last equation at page 74 in iHinzd ( 1975f ) 
shows that the last term in (|3.4p is the global TKE dissipation, 



duj 
dxj 



duj 
dxi 



+ 



dxi 



Equation (|3.4p may therefore be written as 

-Vuv -7 





dui 


dut 




dxj 




- a 













(3.5) 



(3.6) 



The balance in p.6p simply states that the global TKE engendered by the production 
terms, Puv and Vvw, is dissipated into heat by the turbulent viscous stresses. Figure S] 
shows the wall-normal profiles of the three terms whose integrals compose the balance 
p.6p . It is observed that the integrand of Vuv is suppressed near the oscillating wall and 

its peak moves upward. These changes are attributed to the increase of d (^U^ /dy in the 
outer region and to the near- wall reduction of (uv) , as shown in figures [T] (left) and [5] 
(left), respectively. In contrast, Vj- decreases near the wall, but it significantly increases 
at y w 10. As shown in ^ this is directly linked to the enstrophy production through 
the stretching of vorticity fluctuations by the Stokes layer. 



3.3. Total kinetic energy balance 

By summing the global transport equations for MKE, 
balance for the total mechanical energy is found 



w 



and TKE, ([SH), the global 



(3.7) 



The energy input UhT^ (per unit area and unit time), which drives the flow along x, and 
the energy S^,, spent to enforce the wall motion, are dissipated into heat through the 
viscous action of the mean streamwise and spanwise flow gradients, denoted by 1?^/ and 
TDw respectively, and through the visc ous dissipation Vj- of the turbulent fluctuations. 
Note that, as shown bv lLaadhari ( 2007 ) for the uncontrolled flow, Vj- ^ Vjj as Rct oo. 

Figure [5] summarizes and quantifies the global energy balance. The two boxes represent 
MKE and TKE; MKE-x and MKE-z indicate the portion of the MKE balance pertaining 
to the streamwise and spanwise directions, respectively. The light grey portions of arrows 
indicate the energy terms in the fixed-wall case, while the dark grey arrows or portions 
of arrows denote the energy transfers due to the wall motion. The schematic graphically 
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Figure 4. Wall- normal profiles of the integrands oiVuv, Vvw, and Vt in (|3.6p . Time-averaged 
values of production term 3 in H3.3|) . denoted by triangles, of production term 4 in H3.3|) . denoted 
by squares, and of pseudo-dissipation term 6 in p.3p . denoted by circles, for fixed- wall (thin 
lines) and oscillating-wall (thick lines) cases. 




Figure 5. Schematic of the global energy balance for the total mechanical energy. The numbers 
indicate the magnitude of the terms and the additional contribution due to the wall oscillation. 
The light grey portions of arrows denote the contributions at fixed-wall conditions, while the 
dark grey arrows or portions of arrows indicate the changes due to the wall motion. 

highlights that the production terms Vuv and Vvw only transfer energy "internally" 
between MKE and TKE, therefore disappearing from the total energy balance (|3.7p . 

As C/ and Ub are related, the aim is to study how the wall motion acts on Ub to discern 
information on drag reduction. The total energy balance (|3.7p is therefore analyzed in 
more detail because it contains Ub explicitly. As shown in figure El it is first noted that 
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the terms in p.2p pertaining to the streamwise and spanwise directions are decoupled, 
so that 



(3.8) 



and, correspondingly, £w = Vvw + 'Dw- The two terms containing W , i.e. £w and Vw, 
almost balance each other; the difference, given by Vvw (which is absent in (|3.7p ). is 
much smaller than the other terms in (j3.7|) . To gain insight into the changes of Ub, one is 
thus led to investigate how the wall oscillation affects the dynamics of the two remaining 
relevant terms, i.e. Vu and Vq-. The relative contribution of Vjj to the input power UbTw 
increases in the oscillating-wall case. In the fixed-wall case, the input power in viscous 
units is 15.9, 59% (i.e. 9.4/15.9) of which is dissipated by Vu. When the wall oscillates, 
this share in creases to 62% (i.e. 12.1/19.4). Thi s fact agrees with previous studies on 
flow control ( BewlevI l2009l : iFukagata et al. 2009f) . which show that, as C/ decreases as 
the flow tends to the laminar regime, the input power is dissipated more by Vu and less 
hyVr. 

Three different scenarios might explain how the wall motion acts on Pfj and why the 
relative contribution of Vij in the global balance increases during the wall motion, 
(a) In the first scenario, the mean spanwise shear may work directly on Vij. The 

transport equation for 
reads 



(j)U/dyj , the integrand of Vjj (see p.2p ). is thus studied. It 



ld_ 
297 




d^uv dU 
dy^ dy 



1^ 
2d^ 




9y2 • 



(3.9) 



The spanwise velocity W does not appear in the transport equation for (dU/dy)'^, which 
proves that the oscillating wall does not influence the dynamics of the mean streamwise 
flow directly and that the increase of Vu is linked to the modification of turbulent 
dynamics. 

(b) In the second scenario, the spanwise viscous effects directly damp Vuv, so that the 
relative contribution ofVu in (|3.8|) is larger. As Vuv depends on uv and dU /dy (see <\?>.2\ ) 
and the direct action of W on dU /dy has been excluded, the focus is on the transport 
equation for the Reynolds stresses uv (which also appear in the r.h.s. of 



d {uv) 
dr 



_'7l2 



dy 



dy 



dp dp 
dx dy 



^uv du dv 



dy'^ 



dxj dxj 



The mean flow W does not appear explicitly in p.lOp . which demonstrates that the 
oscillation does not work directly on uv either. Therefore, this scenario is excluded as 
Vuv is not immediately affected by the large-scale spanwise flow because neither v/u nor 
dU /dy are. 

(c) The third scenario is rather counterintuitive. The wall oscillation enhances the 
turbulent dissipation, so that the turbulent activity drops due to the increased dissipative 
nature of the flow. A relatively lower level of dissipation of turbulent energy into heat 
thus sets in to balance this new condition. A similar behaviour has been observed in 



applying the suboptimal control theory (|Lee et aLlll998[ ) when the turbulent dissipation 
is employed as a cost function. 

Elucidating the third scenario, i.e. understanding how the wall motion affects the 
dynamics of the turbulent viscous dissipation I?7-, then becomes the aim of the next 
section. 
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4. Turbulent enstrophy balance 

By use of the fluctuating vorticity, term 6 in (|3.3p becomes (|PoDdl2000l) 



dxj dxj * * dxidxj 



The global dissipation of TKE in (13. 3|) becomes 

Vr = , (4.2) 

which follows from the substitution of (|4.ip into (|3.5p . from the periodicity along the 
homogeneous x and z directions, from the velocity fluctuations being zero sd y — and 
because dv'^/dy = at y = h. Note that the viscous dissipation of the total mechanical 
energ y, Vri+Vw+Vq- , equals the global enstrophy only in the case of stationary bound- 
aries ( Davidso"nll2004( l. and therefore not in the oscillating- wall case. However, equation 



(|4.2p is valid for the wall-oscillation case because the turbulent fluctuations vanish at the 
walls. 

Instead of considering the transport equation for the turbulent energy dissipation, we 
opt to study the turbulent enstrophy equation. In the second part of the Appendix, it is 
shown that the form of the turbulent dissipation equation is similar to that of the enstro- 
phy equation and that the dominant terms brought about by the wall oscillation in these 
equations have the same order of magnitude. Expressing Vj- in terms of the turbulent 
enstrophy is more compact than if the turbulent dissipation is used (compare (|4.2I) with 
p.Sp ). The enstrophy equation has the further advantage over the dissipation equation 
that the turbulent pressure does not need to be computed. (This advantage is shared by 
the Orr-Sommerfeld and vorticity formulations of the Navier-Stokes equations over the 
framework involving primitive variables.) Moreover, the physical meaning conveyed by 
the enstrophy equation is arguably more immediate than the one provided by the dissi- 
pation equation; for example, terms 2 and 3 in the turbulent enstrophy equation (j4.3p 
denote production of vorticity, while the corresponding terms in the dissipation equation 
indicate production of turbulent dissipation. 

4.1. Balance equation for the turbulent enstrophy 



The transport equation for the turbulent enstrophy (|Tennekes fc Lumlevlll972l ) reads 




- ^'^^^ + ""^ - -2d-y (^"'''^') (4.3) 

6 7 8 ' 9 

1 d'^ujiuJi duii dwi 

2 dy"^ dxj dxj 
10 n ' 

Term 1 indicates the time rate of change of the turbulent enstrophy. Terms 2 and 3 are 
the production (or removal) of turbulent vorticity caused by stretching (or squeezing) of 
vorticity fluctuations by the mean flow gradients dU /dy and dW/dy, respectively. Terms 
4 and 5 indicate the production of mean and turbulent enstrophy by the stretching of 
fluctuating vorticity through the fluctuating strain rates du/dxj and dw/dxj, respec- 
tively. Terms 6 and 7 represent the exchange of fluctuating vorticity between the mean 
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and the turbulent enstrophy due to the gradients of strcamwise and spanwise mean vor- 
ticity, respectively. They are analogous to the turbulent kinetic energy production terms 
in the MKE and TKE equations p.ip and (13. 3p . Term 8 is the production of turbulent 
enstrophy by stretching of turbulent vorticity through turbulent velocity gradients. Term 
9 denotes the transport of turbulent enstrophy by the fluctuating wall-normal velocity 
component. Term 10 is the viscous transport of turbulent enstrophy and term 11 is the 
viscous dissipation of turbulent enstrophy. The only terms in (|4.3p that become null when 
(|4.3p is made global are term 1, when time averaged because of time periodicity, and term 
9 when integrated along y. 

In contrast to the case of the transport equation (|3.9p for (^dU /dyj , which contributes 

to Vu (see (I3.2p ). W appears exphcitly in terms 3, 4 and 6 of (|4.3p . These terms arise 
only when the wall oscillates. This indicates that the spanwise motion acts directly on 
the turbulent enstrophy and therefore on the global turbulent dissipation "Df. As this 
quantity increases during the wall motion, it is worth studying how these oscillating-wall 
terms contribute to modify the enstrophy balance and, in turn, Uh through p.7p and Cf 
through dMl). 

Figure [S] shows the profiles of the time-averaged terms in the turbulent enstrophy 
balance (|4.3p for the fixed-wall (top) and oscillating-wall (bottom) cases once the new 
fully-developed regime has established. The numbers refer to the terms in (|4.3p and 
the thick lines in the bottom graph highlight terms only occurring duri ng the wall mo- 



tion. The fixed- wall profiles sho w very good agre ement with the ones in lAntonia fc Kim 
( 19941) .iGorski et all ([1994) and I Abe et all (|2009^. (Note that in lGorski et oLI (| 19941 ) a^ 



Abe et al\ (|2009l ) the terms are multiplied by a factor of 2.) In the oscillating-wall case, 
the vorticity production term 3, (cJ^ydW / dyj , is dominant in the proximity of the 
wall, y < 10, over terms 4 and 6, and over the production and transport terms already 
present in the fixed-wall case, i.e. terms 2, 5, 7, 8, 10. This is the key term producing 
turbulent enstrophy (and dissipation). Its physical meaning is further addressed in ij4.3l 
It peaks at y ~ 6 and distinctly affects term 11, the dissipation of turbulent enstrophy, 
at the edge of the viscous sublayer and in the lower part of the buffer region, as clear 
from the similar shapes of the profiles for 2 < y < 20. 

In a very thin near-wall layer, y < 2, term 3 is small. Term 10, the viscous transport 
of turbulent enstrophy, is instead responsible for the intense increase of dissipation of 
turbulent enstrophy there. While the production term (uJ^ydW / dy^ emphasizes the 
direct action of the spanwise shear layer on the turbulent enstrophy, the increase of 
the production term 2, (uj^ydU / dy^ , outlines the indirect effect of the wall motion 

caused by the increment of dU /dy. Term 3 is primarily dominant near the wall, whereas 
term 2 increases at higher locations. This is because term 3 is dictated by the near- 
wall spanwise velocity W , while dU /dy only varies significantly for y > 15, the wall- 
shear stress being constant (see figure [1] (left)). We finally note that the production 

term 5, (ujjdw/dxjdU /dy'^, decreases substantially in the oscillating-wall case, while the 
production term 8 and the transport term 9, which only involve fluctuating quantities, 
are largely unaffected. 

4.2. Transient response of turbulent enstrophy 

Although it is clear from the discussed results that the spanwise shear layer W enhances 
the turbulent dissipation, it is difficult to verify the third scenario proposed in i j3.3l 
through the energy balance shown in figure O i.e. that the relative contribution of Vj- to 
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the global balance decreases. Furthermore, by studying the global balance no information 
is gained on why the TKE decreases when the walls move. We therefore stud y the effect of 
the os cillation at very short times, in line with lQuadrio fc Riccol ( 2003 ) and Xu fc Huane 
(l2005l) . who have used the same approach to study the Reynolds stress budget. 

The temporal evolution of the space-averaged, y-integrated turbulent kinetic energy, 
turbulent enstrophy, squares of the turbulent vorticity components, and terms 2 and 3 in 
(|4.3p are shown in figure[71 where the wall is in motion for t > 0. On this short time scale. 
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Figure 7. Temporal evolution of the space-averaged, y- integrated turbulent kinetic energy (thin 
dashed line), turbulent enstrophy (thin solid line), three vorticity terms forming the turbulent 
enstrophy {ljx denoted by squares, Uy denoted by diamonds, and cj^ denoted by triangles), and 
term 2 (thick dashed line) and term 3 (thick solid line) in (|4.3|) . The oscillation starts at t = 0. 
Terms 2 and 3 are multiplied by a scale factor of 30, while the kinetic energy is divided by a 
scale factor of 100. 




Figure 8. Temporal evolution of the space-averaged, t/- integrated terms in the turbulent 
enstrophy equation (|4.3I) . Thick lines highlight terms only occurring during the wall motion. 
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terms 2 and 3 show an oscillating behaviour, whose period matches well that of the wall 
forcing. Upon the beginning of the oscillation, term 3, denoted by the thick solid line, 
grows abruptly until t — 25, i.e. at a quarter of the oscillation period. It gives a transient 
production of turbulent enstrophy, and, specifically, ofwz, whose production is related 
to the spanwise tilting of Wj, by the Stokes layer, ujydW /dy, in the transport equation 
of oj^. Term 2, which is indicated by the thick dashed line and is non-zero before the 
start of the oscillations, drops substantially up to t—3Q and then becomes larger than 
term 3. As dU /dy is unchanged on such a short time scale, the decrease of term 2 is 
directly linked to uj^y. This can be interpreted as the wall oscillation causing a change 
of the phase relationship between the low-speed streaks (related to ujy) and the quasi- 
streamwise vortices (related to lux)- This scenario is in l ine with early suggestio ns on 
the effects of the forcing on near- wall coherent structures ( Baron fc Quadriol 19961 ). This 



mechanism is indirect, since W does not appear explicitly in term 2. Term 2 also appears 
in the transport equation for and represents the tilting of ujy along x due to dU /dy. 
Figure [7] indeed shows that the first two instantaneous peaks of term 2 agree fairly well 
with those of the transient evolution of tj^. The long-term behaviour of term 2, namely 
the slight increase and the outward shift of its maximum, is due to the increase of dU /dy 
in the bulk of the channel, caused by the increase in mass flow rate. Figure [5] shows 
the temporal evolution of all the space- averaged, y-integrated terms in the enstrophy 
equation ()4.3p . Term 3 is the cause of the short-term transient changes in the enstrophy 
balance, which is evinced by the resemblance of its temporal history with the one of term 
11, the turbulent enstrophy dissipation, up to t=50. 

The turbulent dissipation is therefore enhanced, which causes the monotonic decrease 
of TKE. This feeds back onto the turbulent vorticity and onto term 3, which are both 
diminished because of the weakened turbulent activity. As a direct consequence of the 
attenuation of TKE, the streamwise mean flow accelerates, thereby increasing the stream- 
wise mean velocity. This is evident from the mean streamwise momentum equation, 

„ dU d'^U duv 

because —11 is constant (and positive) and the convective term dU/dt must be positive 
to counteract the decay of the convective transport due to the Reynolds stresses (which 
are larger than the mean viscous terms d^U/dy^). The TKE continuously decreases 
because, although the turbulent dissipation and production are both attenuated, the 
latter is proportionally smaller. As the streamwise flow accelerates, all the quantities 
decrease up to t « 400. This mirrors the trai isient behaviour of the turbulence under 
constant mass flow rate conditions, studied bv lQuadrio fc Riccol (|2003h : on such a short 



time scale, there is no difference between the two constraints. This is further supported 
by the initial attenuation of the wall-shear stress, dU /dy\y=o, which is an immediate 
consequence of the acceleration of the mass flow rate. This is shown by integrating (|4.4p 
along y, 

y=o 

As —11 is positive and the flow-rate term on the r.h.s. is positive, the wall-shear stress 
must be smaller than its steady-state value during the transient evolution. The value 
of the wall stress eventually re-establishes itself in the new fully-developed regime to 
the value imposed by the constant 11. This is evident in figure El where the long-time 
evolution of TKE and of the mass flow rate are also shown. 
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Figure 9. Temporal evolution of the space-averaged wall friction (dashed line), the TKE (dotted 
line), and the mass flow rate (solid line). The TKE follows the temporal behaviour of the friction, 
whereas it is graphically evident how the unbalanced friction causes a flow acceleration with a 
consequent increase of the mass flow rate. All values are normalized by the average uncontrolled 
case. 



The temporal evolution has helped to clarify the action of the wall motion on the 
turbulence dynamics. In particular, as the oscillation initially enhances the turbulent 
dissipation, it is shown that the turbulent activity is suppressed due to the dissipative 
nature of the flow. In the new quasi-equilibrium regime reached after the long transient 
has elapsed, the flow therefore requires a relatively lower level of turbulent dissipation 
because TKE is lower. The third scenario discussed in t i3.3l is therefore at work and term 
3 is key to explaining the lower contribution of Vj- in the global balance. 

We close this section with the schematic in figure [ini The crucial physical processes 
during the temporal flow evolution from the start-up of the wall motion to the new fully- 
developed regime are shown. This last regime is indicated with 'Drag reduction', although 
it should be recalled that, in the present context, the turbulent drag is eventually unvaried 
by design, and that the effect of the oscillations is to increase the mass flow rate. 



4.3. Physical interpretation of (^uj^y^-^ 

Term 3 in the enstrophy equation l|4.3p . (uj^ydW /dy^, has been found to be dominant 

and largely responsible for the change of global turbulent enstrophy and thus for drag 
reduction. It is positive and therefore indicates a production of turbulent vorticity, i.e. 
the mean spanwise flow shear dW/dy acts on the turbulence structures represented by 
the term cJ^y to increase the turbulent enstrophy. We proceed to investigate the physical 
meaning of such interaction in more detail. 
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Figure 10. Schematic of the physical mechanism leading to skin-hiction drag reduction by wall 
oscillations, as discussed in H4.2I The vertical arrows indicate whether the quantities increase or 
decrease. 
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Figure 11. Wall- normal profile of uJ^y at the phase at which it reaches its maximum (thick 
line). Thin lines indicate the sub-terms in (|4.5p . 



The quantity uj^iOy can be expanded as: 



dv du 
dx dz 
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dv dw 
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du du 
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3c 



du dw 
dy dx 

3d 



(4.5) 



Figure [TT] shows that 3c is the largest contributor to uJ^y, and that the next largest 
sub-term in magnitude is term Sd. This confirms the order-of-magnitude analysis in the 
Appendix. Terms du/dz and du/dy, contributing to term 3c, may be linked separately 
to the dynamics of the turbulent low-speed streaks. In the near- wall region, du/dz marks 
the lateral flanks of the low-speed streaks, i.e. the interfaces of the low- velocity and high- 
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Figure 12. Instantaneous contour plots of it and ujyUJ^, for the fixed-wall case, in the x — z 
plane at j/ = 6. For the sake of clarity, only a fraction of the computational domain is shown. 
Gray shading represents the level of streamwise velocity fluctuations, with white corresponding 
to the maximum values and black to the minimum values; max|u| = 8. Contour lines represent 
values of the quantity ujyUj^; contour levels start from ±0.125, are spaced by 0.25, and negative 
values are dashed. 





Figure 13. Instantaneous contour plots of u and ujyUJz, for the oscillating- wall case, in the x — z 
plane at i/ = 6. From left to right and top to bottom: (^=0, 7r/4, 7r/2, 37r/4. Legend is the same 
of figure [II 



velocity regions, while du/dy is related to the eruption of near- wall low-speed fluid to 
higher locations and to the sweep-like motion of high-speed fluid toward the wall. It is 
also noted that the peak location of 3c matches well that of Vj- in figure SI suggesting 

that the enhancement of Vj- is connected to term (uJ^ydW / di/j . 

Figure [m shows contour plots of streamwise velocity fluctuations and ujyUJz in the x — z 
plane at y = 6 for the fixed-wall case. Low- and high-speed streaks show the characteristic 
streamwise-stretched shape. Low-speed streaks are longer and thinner than the high- 
speed ones. Regions of high magnitude of uJyUJz are concentrated near the wall. At y = 6, 
they appear sporadically and always occur at the sides of high-speed regions. Figure \W\ 
shows the contours plots for the oscillating-wall case at four dif ferent phases, where th e 
characteristic cyclic tilting of the near- wall structures is evident (jQuadrio fc Riccdl2003f ). 

The streaks are less energetic, which confirms the attenuation of ^""^^i shown in figure 
[21 Regions of high Iw^Wzj show an analogous tilting, owing to their relationship with the 
velocity streaks. The number, the amplitude, and the spatial size of the tOyUiz pockets 




strongly increase during the wall motion, in line with the observed intensified enstrophy 
fluctuations. 

The interaction between the large-scale oscillating shear layer and the underlining 
vortical structures can be modeled by analogy with the rapid distor tion theory prob lem 
of a large eddy stretching a smaller blob of vorticity, presented by Davidson! ( 2004 ) at 
page 213. We consider small-scale vorticity structures being stretched and c ompressed 
by the large-scale action of the spanwise layer and we extend the model by IPavidson 
(200J) to include the viscous dissipation effects. It is assumed that the forcing induced 
by the wall motion is more energetic than the turbulent fluctuations and it operates on 
a longer time scale. To simplify the problem further, the focus is on the dynamics at a 
wall-normal location ?/ w 6, where, as shown in figure El the vorticity production given 
by the spanwise layer, i.e. term 3 in (14. 3p . is largely balanced by the viscous enstrophy 
dissipation, i.e. term 11 in (|4.3p . The gradient of the spanwise layer, indicated by G, is 
taken as constant in the small region considered. 

In the y — z plane, the enstrophy dynamics may thus be distilled into the following 
simplified equation 



ld_ 
2dt 



G 



dU!y 

dy 



dujz 
dy 



(4.6) 



where only the terms involving wall-normal gradients are retained amongst the dissipa- 
tion terms because they are dominant as revealed by the order-of-magnitude analysis in 
the Appendix. The terms on the r.h.s. of (j4.6p may be written in matrix form as follows 
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As shown in figure 1141 a set of perpendicular axes 
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x,Xs) may be considered where Xn 



is orientated along the vorticity vector in the y — z plane, uiyz = [uJy] uJz]- The angle a is 
defined between Uy^ and the wall-normal axis y. In the new set of coordinates, the terms 
on the r.h.s. of (|4.t)P are written as 

- „ — [^n]yJ\ o o n — •^nni^rn 
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where = ujy +uj1, Snn = sinacosaG, and Dnn — 1, Dns — 0. The other components 
of the strain-rate tensor, S'ns, Ssn, only contribute to change the direction of the vorticity 
vector, not its magnitude. Equation (j4.6p is written as 

''''^-S.^.l-i'-^y. (4.7) 



2 dt " \dy 



Equation (|4.7p may be integrated by Charpit's method (jGarabedianl 1 1 964f) to give 



w„ = a;„,oef^2C^^e;f_Vf5, (4.8) 

stretching dissipation 



where uin,o is the initial magnitude. The constant 



dujn/dt Xy 



dujn/dy At' 

where Ay represents the dissipative scale along the wall-normal direction, and At indicates 
the time scale of the turbulent fluctuations. 

Equation (|4.8|) shows that the spanwise layer may stretch or compress the turbulent 
fluctuations depending on the sign of its gradient and the orientation of the vorticity 
vector. The spanwise layer works by stretching when sin a cos aG > 0. Its action is null 
when the vorticity vector is parallel or perpendicular to the wall, and maximum when 
either i) G is at its negative peak at that y location and simultaneously Uyz is oriented 
at 7r/4 with respect to the axes and along the first or third quadrant in figure [Til or ii) 
G is at its positive peak and simultaneously u)yz is orientated at 7r/4 with respect to 
the axes and along the second or fourth quadrant in figure [T31 It is further noted that 
the temporal rate of growth or decay of enstrophy is never larger than G, and that the 
exponential attenuation through the viscous effects is more intense in time than space. 

We conclude that turbulent vorticity is produced when sin a cos aG > /3^, i.e. when the 
shear-layer production, ruled by the intensity of the large-scale spanwise shear layer and 
by the orientation of the turbulent vorticity vector, overcomes the viscous dissipation, 
whose dynamics is linked to the time scale and the wall-normal spatial scale of the 
fluctuating vorticity. 

4.4. Drag reduction and production of turbulent enstrophy 

The importance of the enstrophy production term 3 has been revealed. However, our 
conclusions are only based on one flow condition, A = 12 and T = 100, yielding R = 0.31. 
It remains to show whether this result can be generalized to other periods of oscillation. 
To this purpose, additional simulations have been carried out by changing T and leaving 
all the other parameters unvaried. Figure [13 shows that R is linearly proportional to 
uj^ydW/dy up to the conditions where R w 0.30, which corresponds to T = 42. 

The linearity between R and the global value of enstrophy production begins to lose its 
validity at T = Topt- 

The global value of the enstrophy term is shown because R is linked via (|2.2p to the 
change oi Ub, a global quantity as defined in ()2.ip . It is known that when T is larger 



than the characteristic life time of the near-wall turbulent structures (jQuadrio fc Ricco 



2004^ . drag reduction drops as the forcing is slow enough to decouple from the near- 
wall turbulence dynamics. At high T, the oscillating wall is not expected to induce drag 
reduction because the near-wall structures are too slowly affected and thus tend to re- 
establish their natural dynamics between sweeps of the Stokes layer. Indeed, the case at 
T = 500 (not shown in the figure) gives a small R — 0.06, and the linear relationship with 
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Figure 15. Drag reduction as function of 
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at different oscillations periods. For 



small values of T the two quantities are linearly related. 



the global term 3 is lost. On the other extreme, the wall motion becomes ineffective and 
produces small R at small T owing to the limited wall-normal extent of the oscillating 
Stokes layer. 

4.5. A note on the scaling parameter by lQuadrio & Riccd t2004) 

In previous works, attempts have been mad e to identify a forc ing parameter which scales 
linearly with drag reduction. In particular, Choi et ali ( 20021 ) have introduced a param- 
eter that has been shown by iQuadrio fc Ricc o (200^ to relate linearly to the amount 
of skin-friction reduction for periods of oscillations smaller than or comparable with the 
optimum. (At larger periods, the above-mentioned decoupling between the forcing and 
the near-wall turbulence takes place.) This scaling parameter was constructed by com- 
bining a characteristic length scale related to the wall-normal distance at which the wall 
motion affects the turbulent structures and the maximum spanwise acceleration of the 
Stokes layer. It reads: 



S 



2./^ In 



A 
Ath 



exp 




where Ath ~ 1.2 is a threshold velocity and y « 6.3 is a wall-normal distance represen- 
tative of the diffusion of the Stokes-layer viscous effects. 

It is then natural to inquire whether there exists a relationship between S and the 
global enstrophy production term (which also relates linearly to R, as shown in figure 
I15p . here expressed as a new parameter, 



Sn=SniT,A) = 



LOzUJy- 



dW 
dy 



(4.9) 
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The first observation is tliat S can be written as 

S=jln(^^^n,^{y), (4.10) 

where flxm is the maximum streamwise vorticity of the Stokes layer aX y ^y. Relation 
(|4.10p endows S with a more direct and physically relevant meaning as it simply states 
that the drag reduction is linearly proportional to the maximum spanwise shear induced 
by the Stokes layer at constant A, and that such shear is most effective when at work 
at y « 6.3. The fact that S relates well with drag reduction is not surprising in view 
of the scaling analysis based on the enstrophy production. The parameter S in (14.10p 
can be seen as a simplified version of 5„ in (|4.9|) . Although 5„ is more elaborate as it 
possesses a precise physical meaning and involves averaging and wall-normal integration, 
the spanwise shear plays a key role in both expressions. The other point of note is 
that (uj^ydW / dy'^ reaches its maximum at y ~ 6.5 for optimum conditions of drag 
reduction, i.e. almost at the same distance at which the correlation between R and S in 
(j4.1L)p is maximum. 

It is noted that other semi-empirical formulas linking R and the wall oscillation pa- 
rameters have been put forward. Bandyopadhyay ( 2006i) developed a formula also based 
on the effect of the Stokes layer on the near-wall turbulence. The central idea is that the 
Stokes layer cyclically re-orients the near- wall vorticity and the drag reduction is linearly 
related to the sine of the maximum angle of vorticity re-alignment with respect to the 
streamwise direction. This approach can be said to belong to the same family of ours, 
although the analysis in the present paper is based on the modification of the turbulent 
flow statistics (specifically related to the dissipation), while Bandyopadhyay 's physical 
model is more directly inspired by the instantaneous action of the spanwise layer on the 
vortical coherent structures. 



5. Summary 

We have described via a DNS study how harmonic wall oscillations are capable of 
increasing the mass flow rate in a turbulent plane channel flow driven by a constant 
pressure gradient. The uniquely defined inner scaling brought about by the constant 
pressure gradient is exploited to ascertain how the oscillations modify the turbulence 
statistics. By looking at the energy fluxes in global form, it has emerged that the energy 
spent to drive the wall motion almost coincides with the viscous dissipation due to the 
oscillating spanwise layer, the difference taking the form of a small turbulence kinetic 
energy production term. The energy balance shows that the enhancement of energy 
intake due to the increased flow rate is mainly balanced by the combined increment 
of dissipation associated with the mean streamwise velocity profile and the turbulent 
dissipation. It is also revealed that the relative contribution of the latter to the total 
dissipation becomes smaller when the wall moves. 

The spanwise oscillating layer shows a direct effect on the turbulent dissipation, which 
is conveniently expressed as the volume integral of the turbulent enstrophy. The amount 
of drag reduction relates linearly with the volume integral of an enstrophy production 
term induced by the spanwise shear layer. The study of the turbulent enstrophy trans- 
port equation reveals that this dominating enstrophy production term synthesizes the 
stretching of the vorticity lines by the oscillating layer, and therefore enhances the turbu- 
lent dissipation. The analysis of the short-term evolution of the flow after the beginning 
of the wall motion shows that the dissipative nature of the near-wall field is responsible 
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for the attenuation of the turbulence intensity and thus for the increment of the bulk 
velocity. 

The study of the turbulent enstrophy in fully-developed conditions to evince which 
term dominates the physics and the analysis of the flow temporal evolution to discern 
how the new regime ensues can be both useful to investigate other turbulent flows mod- 
ified by external agents, such as boundary layers affected by large-scale Lorentz or Cori- 
olis forces, by wall transpiration, or by large temperature gradients. Furthermore, it 
would be of interest to use the approach based on the t u rbule nt enstrophy to investi- 
gate the traveling- wave flow proposed by iQuadrio et al. ( 20091 ) in drag-reduction and 
drag-increase conditions. 



Appendix. Order-of-magnitude analysis on turbulent enstrophy and 
dissipation equations 

The order of magnitude of the terms arising because of the wall motion in the turbu- 
lent enstrophy equation (|4.3p and in the turbulent d issipation equation can be e stimated 
through an analysis similar to the one carried out bv lTennekes fc Lumlev (1972) at pages 
89 and 90. T wo symbols are adopted, fol lowing the introductory discussion on the use 
of symbols in iTennekes fc Lumlev ( 1972[ ). The symbol ~ denotes a crude approxima- 
tion; it highlights the dependence of the term under scrutiny on the characteristic length 
and velocity scales of the turbulent motion. Upon decomposing an enstrophy term into 
sub-terms containing the fluctuating velocit y components, the symbol O denotes its mag- 
nitude in terms of the dominant sub-term. In lTennekes fc Lumlev! (1972), a generic length 
scale is assumed to describe the mean flow motion and the Taylor microscale is taken 
as the reference length scale for the turbulent fluctuations in all directions, suggesting 
that such an analysis is useful primarily for homogeneous isotropic turbulence. However, 
our interest is on the wall-bounded turbulence dynamics with wall oscillations, which is 
strongly anisotropic. It is therefore necessary to distinguish different length and velocity 
scales along the three Cartesian directions. 

The near-wall turbulent dynamics is characterized by three distinct length scales. The 
length scale of the disturbance along z can be taken as A^; = 0(100), na mely the char- 
acteri s tic sp acing of the low-speed streaky structures ( Kline et al. 1967[ ). As shown by 
Riccol (2004), the streaks spacing increases by about a fifth when R k, 0.3, so that the 
order-of-magnitude estimate is still valid. The streaks length, = 0(1000) for fixed - 



wall conditio ns, is representative of the disturbance flow along x (|Kline et al\ 119671 ) 



Riccol ( 20041 ) has shown that A^, decreases by about a third when R « 0.3. The order 
of magnitude of A^; = 0(1000) is therefore applicable under wall-oscillation conditions. 
The length scale along y for the mean flow is the spanwise boundary layer thickness 5, 
defined here as the distance from the wall where the maximum W equals exp(— 1)A. As 
amply verified (see IChoi et al\ (j2002l ) , amongst many) , W agrees well with the laminar 
solution of t he second Stoke s problem for the flow induced by wall oscillations beneath 
a still fluid ( Batchelor 19671 ). so that the spanwise boundary layer thickness can be ap- 
proximated well by 5 = a/T/tt. For T = 100, S « 5.7, so that it can be assumed that 
d = 0(10). The boundary layer thickness 6 can be taken as the characteristic length 
scale for the near-wall disturbance flow because the oscillating boundary layer affects 
the turbulence in a region close to the wall whose width is comparable with S. This is 
shown in figure [5] by the (uv) profile being markedly affected only for y < 25 and by 
the vw profile reaching its maximum at y ~ 15. The mean- flow length scale becomes the 
length scale of the fluctuations along the direction of the shear also in other shear-driven 
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phenomena, such as the penetratio n of free-stream t urbulence into the Blasius boundary 
layer to form the Klebanoff modes ( Leib et allll999l ). In that case, the wall- normal scale 
of the fluctuations within the boundary layer is the Blasius boundary layer thickness. 
The characteristic length scales along the Cartesian directions can therefore be taken as 
Ax > Ax > 6. As for the order of magnitude of the velocity components near the wall, as 
outlined by iPopd ( 2000l ) at page 283, both u and w show a linear growth near the wall, 
but the coefficient is larger for the streamwise component. The wall- normal component 
V is smaller than both u and w because it grows quadratically from the wall. The hy- 
pothesis u > w > V can therefore be adopted. The terminologies 'larger' and 'smaller' 
are used in the order-of-magnitude sense and the time-averaging symbol is omitted for 
brevity. 

The turbulent enstrophy equation (14. 3p is considered first. The order of magnitude of 
terms 3, 4 and 6, arising in (|4.3p because of the wall oscillation, is estimated. Term 3 can 
be first decomposed as follows. 
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and the order of magnitude of each sub-term is 
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It is evident that term 3c, {du/dy){du/dz)dW /dy, is dominant. It follows that 

Term 3 : uj^Uy—— = C -r^— 

It further occurs that term 3d > term 3a > term 36. The magnitude of term 4 is estimated 
as follows. 
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Terms 4a and 4c are larger than 46 and 4c? because w > v and (5 < Az, so that 



Term 4 : LUi 



du dW ^ / uwA\ 
dxi dy \5'^\x) 



(5.2) 



Note that this represents an upper bound because terms 4a and 4c may add to produce 
a term of the order of magnitude given in (|5.2p or give a term of smaller amplitude if 
these terms are of opposite sign. The magnitude of term 6 can be estimated as follows. 



Term 6 
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The term —vdw/dy{d'^W/dy^) is clearly dominant because w > v and 5 > Az. It follows 
that 



Term 6 
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Qy2 y p 

In order to compare term 6 with term 3, we resort to the continuity equation, as follows 

dv du V u 
dy dx S Aj, ' 



Term 6 
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Since u > w and Xx > A^, one obtains 
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> Term 6 : O 



•jA 



uwA 



Terms 4 and 6 are either comparable, when the upper bound case for the order-of- 
magnitude estimate for term 4 is considered, or term 4 < term 6 if the two comparable 
leading terms in 4 have opposite sign. It can be concluded that term 3 > term 6 > term 
4, which is the result found through the numerical simulations. 

The transport equation for the turbulent energy dissipation, called e here 

_ dui f dui duj 
dxj \dxj dxi 

is now studied I'm ansour et a/.lll989l : lFischer et a/.ll2001 ). For the case of turbulent chan- 
nel flow with spanwise wall oscillations, the equation reads 
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The order of magnitude of the terms arising in (|5.3p because of the wall motion can 
be estimated through an analysis similar to one for the enstrophy equation (j4.3p . The 
magnitude of term 3 in (|5.3p is found as follows. 



\ 
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Term 36 is dominant, so that 



Term 3 



dw dv dW 
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The magnitude of term 5 in (|5.3p is estimated as foUows 
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Term 5a is dominant. It follows that 

dui dui dW 



It is found that 
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Term 5a is estimated to be the largest one amongst the terms in (|5.3p induced by the 
wall motion. This result confirms the analysis of the turbulent enstrophy, where term 3, 
of the same order of magnitude, emerges as dominant. 
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